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Resumen 


Se presenta un paquete de cómputo para análisis de frecuencias de caudales 
máximos anuales basado en la aplicación de las hojas de cálculo de uso común 
provistas por Excel” (Excel es una marca registrada de Microsoft Corporation, 
Inc.), el cual usa diversas funciones de distribución de probabilidad. El paquete 
de cómputo maneja la estimación de los parámetros, de los valores de diseño y 
de sus límites de confianza para varios periodos de retorno por diversos métodos 
de estimación de parámetros, valores de diseño y límites de confianza. El paquete 
de cómputo usa hasta tres medidas de bondad de ajuste y produce varias gráficas 
que ayudan a la comparación visual entre los modelos en competencia. Se incluye 
en el artículo un ejemplo de aplicación, utilizando la distribución de probabilidad 
general de valores extremos, y el método de estimación es el de momentos de 
probabilidad pesada, para demostrar la facilidad de uso del paquete de cómputo. 


Palabras clave: análisis de frecuencias de caudales máximos anuales, hojas de 
cálculo, funciones de distribución de probabilidad, estimación de parámetros, 
estimación de valores de diseño, límites de confianza. 


g) Valores Extremos Tipo I (VED. 
h) General de Valores Extremos (GVE). 


Introducción 


En la planeación y el diseño de sistemas 


hídricos, el análisis de frecuencias de caudales 
máximos anuales tiene una importancia central. 
Por las características de los valores de diseño, 
que están ligados a un periodo de retorno, 
es obligado el uso de modelos matemáticos, 
llamados en este caso funciones de distribución 
de probabilidad. Entre las 
distribución de probabilidad más usadas en el 
análisis de frecuencias de caudales máximos 
anuales están las siguientes (Kite, 1988; Rao y 
Hamed, 2000; Raynal-Villaseñor, 2010): 


funciones de 


a) Normal (NOR). 

b) Log-Normal de dos parámetros (LN2). 
c) Log-Normal de tres parámetros (LN3). 
d) Gamma (GAM). 

e) Pearson Tipo III (PITD. 

f)  Log-Pearson Tipo II (LPITD. 


Hay disponibles varios paquetes de cómputo 
para realizar análisis de frecuencias de caudales 
máximos basados en casi todas la distribuciones 
de probabilidad mostradas previamente, como 
lo son el HEC-SSP del U.S. Corps of Engineers, 
US-CE (2010), el propuesto por Reiss y Thomas 
(2007), el propuesto por Rao y Hamed (2000), y 
el propuesto por Kite (1988), entre otros. 

A la luz de las aplicaciones computacionales 
en la educación y entrenamiento en ingeniería 
y ciencia aplicada, se desarrolló un paquete 
de cómputo personal en el área de análisis de 
frecuencias de caudales máximos anuales, con 
énfasis particular en ingeniería hidrológica, pero 
que es fácilmente adaptable a cualquier área que 
use el análisis de frecuencias como herramienta 
estadística de análisis de muestras de datos. El 
paquete de cómputo provee diferentes modelos 
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para ser aplicados y produce una gran cantidad 
de resultados que pueden ser utilizados por el 
usuario, ya sea en forma numérica o gráfica. La 
plataforma computacional escogida fue aquella 
provista por Excel”, lo cual lo diferencia de los 
mencionados previamente y lo hace accesible 
a todos los usuarios de una computadora 
personal. Al programa de cómputo resultante 
se le ha llamado FLODRO 4.0, como será 
referido de aquí en adelante. 

El artículo muestra las principales compo- 
nentes de FLODRO 4.0 e incluye un ejemplo 
de aplicación usando el paquete de cómputo 
personal citado, mostrando los principales 
resultados que dicho paquete de cómputo le 
provee al usuario. 

Las aplicaciones mostradas previamente 
están debidamente documentadas a detalle 
en Raynal-Villaseñor (2010) y algunas de ellas 
pueden consultarse en Kite (1988), y Rao y 
Hamed (2000) En Raynal-Villaseñor (2010) 
se indica cómo se pueden obtener las hojas 
de cálculo de Excel” que se han desarrollado 
dentro de FLODRO 4.0. 


Estructura del paquete de cómputo FLODRO 
4.0 


El paquete de cómputo FLODRO 4.0 ha sido 
desarrollado dentro del ambiente de cómputo 
provisto por Excel? (Excel es una marca 
registrada de Microsoft Corporation, Inc.). El 
modo interactivo en que se usa FLODRO 4.0 
lo hace muy amigable con el usuario, como es 
típico en el ambiente de Windows” (Windows 
es una marca registrada de Microsoft Cor- 
poration, Inc.). 


En cualquier etapa, el usuario tiene todo el 
control sobre los procesos que está ejecutando, 
desde la entrada de datos hasta la impresión 
de resultados. El paquete de cómputo personal 
FLODRO 4.0 tiene la estructura mostrada en la 
figura 1. Todas las funciones de distribución 
de probabilidad citadas en la sección anterior 
están programadas en FLODRO 4.0. 

El paquete de cómputo personal FLODRO 
4.0 puede realizar los siguientes procesos, 
como se muestra en la figura 2: 


a) Estimación de parámetros por los méto- 
dos de momentos, máxima verosimilitud 
y momentos de probabilidad pesada, 
cuando este método está disponible para 
todas las funciones de distribución citadas 
anteriormente. Se produce una gráfica que 
resume estos resultados. 

b) Pruebas de bondad de ajuste, como el 
error estándar de ajuste (EEA) (Kite, 
1988); desviación media absoluta relativa 
(DMAR) (fain y Singh, 1987), y el Criterio 
de información de Akaike (CA), sólo en 
el caso del uso del método de máxima 
verosimilitud (Akaike, 1974). 

c) Cálculo de los eventos de diseño para 
varios periodos de retorno, por los méto- 
dos de momentos, máxima verosimilitud 
y momentos de probabilidad pesada, 
cuando este método está disponible, para 
todas las funciones de distribución citadas 
anteriormente. 

d) Estimación de los límites de confianza para 
los eventos de diseño, por los métodos 
de momentos, máxima verosimilitud 
y momentos de probabilidad pesada, 


Funciones de distribución 
de probabilidad 


FLODRO 4.0 


Figura 1. Estructura del paquete de cómputo FLODRO 4.0. 
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Función de distribución 
de probabilidad 


Comparación entre histograma y modelo 


Figura 2. Opciones de análisis del paquete de cómputo FLODRO 4.0. 


cuando este método está disponible, para 
todas las funciones de distribución citadas 
anteriormente. Se produce una gráfica que 
resume estos resultados. 

e) Una comparación gráfica entre la mejor 
opción de los modelos producidos por 
cada uno de los métodos de estimación de 
parámetros antes citados y el histograma 
de la muestra de datos. 


El paquete de cómputo personal FLODRO 
4.0 ha sido diseñado para tener requerimientos 
mínimos de memoria y de equipos periféricos 
ligados a la computadora personal. 

Las gráficas que produce FLODRO 4.0 
pueden ser impresas en cualquier impresora 
común. En el caso de las gráficas que muestran 
los modelos y los límites de confianza de los 
valores de diseño, se ha usado en el eje de las 
abscisas el valor reducido de Gumbel, para 
producir un papel de Gumbel electrónico, 
dadas las limitaciones que tiene Excel? en 
cuanto a poder manejar un eje con escala 
probabilística o con periodos de retorno, como 
puede observarse en las figuras 7 y 8. 

Estas particularidades hacen del paquete 
de cómputo personal FLODRO 4.0 muy 
apropiado para la educación, entrenamiento y 
educación continua en el área de la ingeniería 


hidrológica en particular, y de la ingeniería y la 
ciencia aplicada en general, principalmente en 
países en vías de desarrollo. 


Ejemplo numérico de aplicación 


La estación hidrométrica Huites está localizada 
en el río Fuerte, en el noroeste de México, y 
ha sido seleccionada para hacer el análisis de 
frecuencias de caudales máximos anuales, 
usando la distribución de general de valores 
extremos (GVE) y el método de momentos de 
probabilidad pesada (MPP), según Raynal- 
Villaseñor (1987), para estimar sus parámetros, 
valores de diseño y límites de confianza. La 
ubicación geográfica de la estación hidrométrica 
Huites, México, se muestra en la figura 3. 

El primer paso es obtener los estadísticos 
más relevantes de la muestra de caudales 
máximos anuales y que algunos de ellos serán 
usados posteriormente, los cuales se muestran 
en la figura 4. 

La función de distribución y de densidad de 
probabilidad para la distribución GVE es como 
sigue (Raynal-Villaseñor, 2010): 


1/8 


(1) 


F(x) = exp+- 1 -E%0)B 
al 
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Figura 4. Datos estadísticos de la muestra proporcionados por FLODRO 4.0. 
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a (2) 


” 


donde F(x) y f(x) son las funciones de dis- 
tribución y de densidad de probabilidad x, 
respectivamente; x, a. y B son los parámetros 
de ubicación, escala y forma de la distribución 
GEV. El dominio de la variable aleatoria es x, + 
a/B<x<oo, 

En la distribución GVE se estiman sus 
parámetros por el método de momentos de 
probabilidad pesada, por medio del método 
propuesto por Raynal-Villaseñor (1987), de la 
siguiente manera (Raynal-Villaseñor, 2010): 


a) Parámetro de ubicación: 


A 
A 


+o-Mp+z[r+8)-1] (3) 


b) Parámetro de escala: 


c) Parámetro de forma: 


Ñ B 
F(B) =2 (Cp +23) -(Cay+3)=0 (5) 


donde: 
Cs Ida abi (6) 
Mo Mi 


Los momentos de probabilidad pesada 
requeridos por las ecuaciones (3)-(6) son como 


sigue: 


- 1 Y 
Mo) da 2 me (7) 


Mau | x,(N -i) (8) 


Las pruebas de bondad de ajuste selec- 
cionadas para la distribución GVE son: 


a) Error estándar de ajuste (EEA) (Kite, 1988): 


Ni 1/2 
S (xi ya? 


EEA= MEN” 2) (10) 


b) Desviación media absoluta relativa 
(DMAR) (ain y Singh, 1987): 


(x,-y;) 


Xx; 


N 
DMAR- LS (10) 
N 


i=l 


Los valores de diseño y sus límites de con- 
fianza son calculados de la siguiente manera: 


P 
oro runa) | (12) 


Q=Q7 +2, Sy (13) 


(14) 
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Las derivadas parciales anteriores tienen la 
siguiente estructura: 


“T-=1 (15) 


B 
ox _1 a] (16) 
B T 


B 
Ma -ta(1-2 
2 all d 


as 
9B B p 
1 1 
+|-En[1-— || |En| -Ln| 1-— 
da L, 


Los coeficientes 10, de la matriz varianza- 
covarianza de la distribución GVE, que apa- 
recen en la ecuación (14), pueden ser evaluados 
a partir de los valores contenidos en el cuadro 
1. 

Al usar las fórmulas anteriores, los pará- 
metros, las pruebas de bondad de ajuste, los 
valores de diseño y sus límites de confianza 
producidos por FLODRO 4.0 son mostrados 
en la figuras 5 y 6. La figura 7 muestra los 
modelos producidos por los métodos de 


momentos, máxima verosimilitud y momentos 
de probabilidad pesada en la estimación de 
parámetros de la distribución GVE, basados 
en la muestra de caudales máximos anuales 
de Huites, México. En la figura 8 se muestra la 
representación gráfica de los valores de diseño 
y sus límites de confianza para el modelo, que 
fue escogido de acuerdo con los indicadores de 
bondad de ajuste manejado en el paquete de 
cómputo. La comparación entre el histograma 
de la muestra de datos de caudales máximos 
anuales y la función de densidad teórica ajus- 
tada se muestra en la figura 9. 


Conclusiones 


Se ha presentado un paquete de cómputo 
personal para el análisis de frecuencias de 
Debido a los 
requerimientos mínimos de memoria y de 


caudales máximos anuales. 


accesorios periféricos, como ha sido demos- 
trado en el artículo, este paquete de cómputo 
personal es muy útil en la educación y 
entrenamiento en ingeniería hidrológica, y 
donde el usuario tenga acceso a una compu- 
tadora personal y una impresora. El uso de 
Excel”, una herramienta computacional de 
uso común, como el ambiente computacional 
seleccionado por el paquete de cómputo, lo 
hace asequible a cualquier usuario con acceso 
a una computadora personal. 


Cuadro 1. Coeficientes 1,, de la matriz varianza-covarianza de los estimadores MPP de los parámetros de la 
distribución GVE (Hosking et al., 1985). 


a 

S B w,, Do w,, w, w,, w,, 
E -0.4 1.6637 1.3355 1.1405 1.8641 1.1628 2.9092 
a -0.3 1.4153 0.8912 0.5640 1.2574 0.4442 1.4090 
E -0.2 1.3322 0.6727 0.3926 1.0013 0.2697 0.9139 
E -0.1 1.2915 0.5104 0.3245 0.8440 0.2240 0.6815 
E 0.0 1.2686 0.3704 0.2992 0.7390 0,2247 0.5633 
3 0.1 12551 0.2411 0.2966 0.6708 0.2447 0.5103 
_ 0.2 1.2474 0.1177 0.3081 0.6330 0.2728 0.5021 
El 0.3 1.2438 -0.0023 0.3297 0.6223 0.3033 0.5294 
$ 0.4 1.2433 -0.1205 0.3592 0.6368 0.3329 0.5880 
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4 
5 
5] 
7 2498.9627 
8 Miy= 719.5557 m/s Cry -1.2770 
A A 
10 
1 
12 Valor de prueba de f = -0.3894 1+P= 0.6106 
13| F(p)= -4.84118-07 T (1 +f)= 1.4654 
14 
15 Parámetro de ubicación 1412.6485 m?/s 
16. Parámetro de escala 908.9530 m?/s 
17 Parámetro de forma -0.3894 
18 
19| EEA= 601 DMAR = 7 
Figura 5. Estimación de parámetros y sus medidas de bondad de ajuste (GVE-MPP) para la estación 
hidrométrica Huites, México. 
[al A B a D E € G H 
21 
2| w,= 1.6269 1, 1.2741 ,3= 1.0538 
23 w,= 1.7630 W,,= 1.0521 y = 1.7497 
24 | var(x)= 26355.9450 cov (xy a)= 20639.9854 
225| var(a)= 28559.8526 cov (xy P) = 18.7808 
26| var(P)= 0.0343 cov (xa, PB) = 0.0206 
27 
28 
29 | Periodo de Raiz cuadrada Límite inferior Valores de |Límite superior 
30 retorno del errorest. de confianza diseño de confianza 
31| Tr S, 95% Q, 95% 
32 (años) FO) = dx, / da dx / 98 (m?/s) (m/s) (m/s) (m/s) 
33 2 0.5 0.3940 -67.1824 29.5992 1713 1771 1829 
34 | la] 0.8 2.0373 -1523.2495 70.1988 3127 3264 3402 
35 10 0.9 3.6003 -4213.2392 139.7205 4411 4685 4959 
36| 20 0.95 5.5960 -8978.6111 264.9121 5980 6499 7018 
87 | 50 0.98 9.1668 -20222.1335 561.3453 8645 9745 10845 
38 | 100 0.99 12.8331 -34441.4956 935.7002 11243 13077 14911 
“39 | 500 0.998 26.2995 -101650.9497 2698.8870 20028 25318 30608 
40 1000 0.999 35.2514 -155159.1684 4099.1379 25420 33455 41489 
41] 5 000 0.9998 68.2202 -388773.9389 10199.9060 43430 63422 83413 
42 10 000 0.9999 90.1555 -565810.3018 14816.9580 54319 83360 112401 
| Límites de confianza de dos colas: La= 1.96 


Figura 6. Eventos de diseño y de sus límites de confianza (GVE-MPP) para la estación hidrométrica Huites, México. 
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Análisis de caudales máximos anuales 
Estación hidrométrica: Huites, México 
(1942-1992) 

Distribución General de Valores Extremos 


Q (m8/s) 


15  -05 0.5 1.5 2:5 3.5 45 55 06.5 
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Figura 7. Curvas de frecuencias teórica y empírica para la estación hidrométrica Huites, México. 


Análisis de caudales máximos anuales 
Estación hidrométrica: Huites, México 
(1942-1992) 


Límites de confianza: distribución GVE (MPP2) 
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Figura 8. Curvas de frecuencias teórica y empírica, y sus límites de confianza para la estación 
hidrométrica Huites, México. 
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Histograma 
Estación hidrométrica: Huites, México (1942-1992) 
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Figura 9. Histograma y la función de densidad de probabilidad teórica para la estación hidrométrica Huites, México. 
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Nomenclatura 

a parámetro de escala. 

B parámetro de forma. 

Sy raíz cuadrada del error estándar. 

o? varianza. 

T. periodo de retorno. 

K, factor de frecuencia. 

x, parámetro de ubicación. 

LV función logarítmica de verosimilitud. 


%. datos de muestra. 

N número de datos de muestra. 

Y, valores producidos por el modelo. 

Q, límite de confianza inferior o superior. 
Za valor normal estándar. 

Q valor de diseño. 
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Abstract 


RAYNAL, J.A. FLODRO 4.0: A personal computer model for flood frequency analysis. Water 
Technology and Sciences (in Spanish). Vol. IL No. 2, April-June, 2012, pp. 41-50. 


A computer package is presented for flood frequency analysis based on the application of 
commonly used Excel spreadsheets (Excel is a registered trademark of Microsoft Corporation, 
Ínc.), using several probability distribution functions. The computer package handles the 
estimation of parameters, design values and confidence limits for several return periods using 
several methods to determine these estimations. The computer model uses up to three goodness 
of fit tests and is able to produce several graphic displays that aid visual comparisons among 
competing models. The paper includes an application example using the generalized extreme 
value probability distribution function and probability weighted moments for the method of 
estimation, to show the ease of using this computer package. 


Keywords: flood frequency analysis, spreadsheets, probability distribution functions, esti- 
mation of parameters, estimation of design values, confidence limits. 
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